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ABSTRACT 

An updated grid of stellar yields for low to intermediate-mass thermally-pulsing 
Asymptotic Giant Branch (AGB) stars are presented. The models cover a range in 
metallicity Z = 0.02, 0.008, 0.004, and 0.0001, and masses between 1M to 6M Q . New 
intermediate-mass (M > 3Af Q ) Z = 0.0001 AGB models are also presented, along 
with a finer mass grid than used in previous studies. The yields are computed using 
an updated reaction rate network that includes the latest NeNa and MgAl proton 
capture rates, with the main result that between ~ 6 to 30 times less Na is produced 
by intermediate-mass models with hot bottom burning. In low-mass AGB models we 
investigate the effect on the production of light elements of including some partial 
mixing of protons into the intershell region during the deepest extent of each third 
dredge-up episode. The protons are captured by the abundant 12 C to form a 13 C 
pocket. The 13 C pocket increases the yields of 19 F, 23 Na, the neutron-rich Mg and Si 
' isotopes, 60 Fe, and 31 P. The increase in 31 P is by factors of ~ 4 to 20, depending on 

\ the metallicity. Any structural changes caused by the addition of the 13 C pocket into 

the He-intershell are ignored. However, the models considered are of low mass and any 
such feedback is likely to be small. Further study is required to test the accuracy of 
the yields from the partial- mixing models. For each mass and metallicity, the yields 
£N| ■ are presented in a tabular form suitable for use in galactic chemical evolution studies 

or for comparison to the composition of planetary nebulae. 
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1 INTRODUCTION ing on the core mass. The energy from the thermal pulse 

_ . . „. _ , . , . (TP) or flash drives a convective pocket in the He-rich in- 

lhe Asymptotic Giant Branch (AGB) phase is the last nu- , nl , , . r TT , . , 

, „ . , tershell, that mixes the products ol He- burning nucleosyn- 

clear burning phase for stars with initial masses between . . . . „ „ . jn 

:f. thesis withm this region, hollowing a IP, the convective 

~ O.oMf^ to 8M^, where the exact limits depend on the , , ,. > . . , 

^ y ~> envelope may move inward (m mass) to regions previously 

initial metallicity, Z. During the AGB there is a complex in- , , „ , , . ,. , , mi ■ ■ j 

. , , ° , r mixed by the flash-driven convective pocket, this inward 

terplay of nucleosynthesis and mixing that alters the surface , , . , . , , , , , . , 

f movement of the convective envefope is known as the third 

composition of the star. The enriched AGB envelope is even- , , . . .. . r ... 

, . , .. , ln „. , , dredge- up ( 1DU), and is responsibfe for enriching the sur- 

tuafly cxpclfcd into the lntcrstcJlar medium (IbM) by a sfow c . \2A j n j r tt u n 

, T „,, . , face m C and other products of He-burning, as weff as 

steflar wind thus enriching the local ISM with the products , , . . , , . . , . ,, . . 

heavy elements produced by the s process m the He-rich 

of hydrogen and helium burning, and heavy elements pro- . , . mr^TT ,i , i ,i 

intershell. hollowing the iDU, the star contracts and the 

duced by the slow neutron capture process (the s process). TT , „ . ., , . n . „ , 

, , , , H-shefl is re-igmted, providing most of the surface fummos- 

Hence these stars are important contributors to the chenu- ., . , . , , . , T a ^t, , .., . ... , 

ity for the next mterpufse period. In AGB stars with initial 

cal evolution of galaxies and stellar systems. For recent re- s .„ ,, , , n 

r- -i masses > ANlp> , the base of the convective envefope becomes 

views of AGB evolution and nucleosynthesis see Busso ct al. , , , . , . /, . 

i 71. ,rrz 1 j u hot enough to sustain proton-capture nucleosynthesis hot 

19991) and lHerwid 20051). , ., ° . mm ^„„ F . * v 

' ~y r ~~' bottom burning, HBB). HBB can change the surface com- 

Bricffy, during the TP- AGB phase the He-burning shefl , , . ,,,,,, 

4 , position because the entire envelope is exposed to the hot 

becomes thermally unstable every 10 years or so, depend- , . . r ,, , , . . , , . , 

J ii- burning region a tew thousand times per mterpufse period. 

AGB stars with HBB have short lifetimes (r < 100 Myr) 

. and are one of the stellar sites proposed as the polluters of 
E-mail: akarakasSSmso.anu.cdu.au 
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globular clusters (|Cottrell fc Da Costalll98ll ; iGratton et al.l 
l2004l ; lRenzinii r2008), even if qua ntitative problems with the 
models exist (e.g.. iFenner et al.ll2004l ). 

In iKarakas fc Lattanziol (l2007f ) we presented AGB 
model data and stellar yields for masses between 1 and 6Mq , 
and for metallicities Z = 0.02, 0.008, 0.004, and Z = 0.0001. 
The main drawbacks of that study were that the reaction 
network dated back to 2003 and there have been signifi- 
cant changes to some of the important proton and a-capture 
rates since that time. Second, the initial compositions for 
the Z — 0.008 and 0.004 models were not scaled solar, but 
had sub-solar C, N, and O compositions that were thought 
appropriate for the Large and Small Magellanic Clouds. We 
also did not investigate the effect of including a 13 C pocket in 
the top of the intershell region, which is required to produce 
neutro ns by the reaction 13 C (a, n) le O during the interpulse 
period l|Straniero et al.lll995l ). The neutrons released are re- 
quired for the s process but can also effect the abundance 
of elements lighter than iron in the He-intershell, and hence 
the stellar yields. Here we include a partially mixed zone 
in low-mass AGB models of 2M , Z = 0.0001, and 3M Q , 
Z = 0.02, 0.008, and 0.004. 

The main aims of this paper a re to provide an up- 
date t o the stellar yields presented in lKarakas fc Lattanziol 
l|2007f ) using newer reaction rates, and to provide yields 
calculated from scaled-solar initial abundances for the 
Z = 0.008 and 0.004 models. We also present yields cov- 
ering a finer mass grid than used previously, and new 
intermediate-mass Z = 0.0001 AGB models. Here we de- 
fine intermediate-mass models as having M > 3Mq, ex- 
cept at Z = 0.0001 where M > 3M . Furthermore, we 
examine the effect of a partial mixing zone on the stel- 
lar yields of low-mass AGB stars. The new results are 
compa red to the yields presented in IKarakas fc Lattanziol 
|2007 ), and to other recent AGB yields in the literature 
(e.gjstancriffe fc Jeffervll2007t IStancliffe fc Glebbeekll200Sl ; 
ICristallo et al.ll2009l ; IVentura fc D'Antonall2009h . Only com- 
parisons between detailed AGB models are made; com par- 
isons to synthetic AG B models (e . g., Marigol l200lh was 
discussed in detail in llzzard et al.l l|2004 ). We also limit 
our discussion to AGB mod els with [Fe/H] > —2.3; see 
ICampbell fc Lattanziol (|2008h for yields from very low- 
metallicity AGB models. 



2 NUMERICAL METHOD 

We calculate the structure first and perform detailed nucle- 
osynthesis calculations afterward, using a post-processing al- 
gorithm. The details of this procedure and the codes used to 
compute the models ha ve been previously de sc ribed in some 
detail , see for example IKarakas et alJ (120021 ) , iLuearo et al.l 
( 2004) and IKarakas fc Lattanziol (^OPTf ). In regards to the 
input physics used in the stellar structure computations, 
we emplo y the mixing-length the ory of convection with 
q = 1.75 ( Frost fc Latta nzio 199 6;). At high temper atures 
we use the OPAL opacities (Iel esias fc Rogers! 1 1996^1 . and 
at low temper atures we employ two differ ent prescriptions. 
The models of IKarakas fc Lattanziol (120071 ) and most of the 
new models here (see below) include an approximate treat- 
ment for the molecular opacities (in particular CN, CO , 



and corrected bv lChiosi etail (|l993h . These fits do include 
some compositional dependence, but do not account for 
large variations in C/O or nitrogen. The second prescrip- 
tion we use is to includ e the low-temperatur e opacity ta- 
bles fr om lFerguson et al.l (|2005l ) in place of the lBessell et al.l 
(|l989h fits. 

The new AGB structure m odels are presented in Sec- 
tion[3] Some of these models have lVassiliadis fc Wood! (|l993l . 
hereafter V W93) mass loss on the AGB , as do the models 
presented in lKarakas fc Lattanziol (|2007T ). We also compute 
new intermediate-mass Z = 0.0001 AGB m odels using the 
Reimer's mass-loss formula (|Reimerslll9"75l . hereafter R75), 
and use them in pla c e of t he VW93 models presented in 
IKarakas fc Lattanziol (|2007h ■ The main reason for this is to 
use a higher mass loss rate than given by VW93 at this 
metallicity, which would allow the 6Mq, Z = 0.0001 model 
to be evolved to ve ry low envelope mass. The VW93 6Mq 
model presented in lKarakas fc Lattanziol l|2007h had ~ 100 
TPs with very little reduction in envelope mass. The VW93 
models are compared to the R75 models in Section [3] 

We scale the R75 mass-loss formula to the metallicity 
of the model according to 



dM 
dt 



V{Z/Z Q ) x Mb 



(1) 



where dM /dt is in Mq yr , Z the global metallicity of the 
stellar model, and Zq the solar metallicity (0.02), and 



M R =* 1O -V* /R0)(L/L0) 



(M/M ) 



(2) 



ment tor tne molecular opacities (m particular PIN, LU , 
H2O, TiO) using the formulations from Be ssell et al.l (|l989l ) 



here R, L and M are the radius, luminosity and mass 
of the star (in solar units), respectively. We note that 
iDecressin et al.l l|2004h apply a similar Z scaling to their 
low-metallicity, intermediate-mass AGB models. The free 
parameter rjn used for each model is provided in Section [3j 
Equation [2] was used on the the first giant branch with 
77R, = 0.4 in all stellar models. 

The new stellar structure models are mostly computed 
with t he same version of the Mona sh stellar structure code 
used in lKarakas fc Lattanziol (2007). The intermediate-mass 
Z = 0.0001 AGB models of 4, 4.5, 5, 5.5, and 6M are com- 
puted us ing an updated versio n of the structure code that 
includes iFerguson et al.l (|2005l ) low-temperature opacities. 
To check for consistency, we compute a 4Mq, Z = 0.0001 
model using both versions of the code. The results of the 
models are presented in Section [31 

The technique used in the post-processing nucleosyn- 
thesis computati ons have been describ ed in detail elsewhere, 
see for example IKarakas etail < |2009h . Here we summarize 
information pertinent to the current models. We assume a 
network of 77 species from hydrogen to sulphur, along with 
a small group of iron-peak elements. Hence this study is lim- 
ited to the nucleosynthesis of light elements. Yields of heavy 
elements produced by the s pro cess using the s ame codes 
are now becoming available, see IKarakas et al.l (120091 ) and 
IChurch et al.l (|2009h . Most of t he 589 reaction rates a re from 
the JINA REACLIB database l|Sakharuk et alJl2006( Fl The 
details of the reaction rates used are described in Section l2~T1 
We assume Z — 0.02 for the solar composition and take the 



http: / /groups. nscl.msu.edu/jina/reaclib/db/ 
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initial abundances from lAnders fc Grevessd (|l989h . For all 
other models we assumed a scaled solar initial composition 
for all species. 

2.1 Updated reaction rates 

The reaction rate network used in the post-processing 
calculations is the rec ommended JINA REACLIB library 
l|Sakharuk et all l2006h with the following changes. We in- 
cluded th e rec o mmen ded rate for t he 18 F(q,p) 21 Ne reac - 
tion from Eel (|2006h and used in iKarakas et aH (|200Sh . 
We used the 22 Nefq, n) 25 Mg and 22 Ne(a,7) 26 Mg rates 
from Karakas et al | l|200rj) . and the 19 F(«,p) 22 Ne rate from 



lllgalde et al l (|200S ) 

Our cur rent reaction rates di f fer si gnificantly to the 
rates used in IKarakas fc Lattanziol (|2007l ) for the following 
reasons. First, the NeNa and MgAl chain proton-capture 
rates have been updated. The mos t important c hange s 
are the 21 Ne(p, 7) 22 Na rate from llliadis et al.l l|200^ . 
the 22 Ne(p, 7) 23 Na rate from lHale et alJ <|2002h, and the 
2 3 Na(p , 7) 24 Mg and 23 Na(p, a) 20 Ne rates from Hale et all 



more C than is left over from CN cycling in the H-shell. 
Hence some mechanism to mix protons from the H-rich en- 
velope into the intershell is needed to produce the extra 13 C. 
In our models, protons are mixed into the intershell region 
by artificially adding a partial mixing zone (PMZ) at the 
deepest extent of each TDU. These protons are captured 
by the abundant 12 C to form 13 C and 14 N, resulting in the 
formation of a 13 C pocket. In the 13 C pocket, neutrons are 
liberate d by the reaction 13 C (a, n) le O during the interpulse 
period (|Straniero et alj[l995f ). 

The timescale for neutron production and the neutron 
source determine the resulting s-process element distribu- 
tion. The details of how the 13 C pocket forms and its extent 
in mass in the He-intershell are still unknown, although var- 
ious mechanisms have been proposed inclu ding convective 
overshoot, rotation, and gravity waves; see Irlerwid (120051 ) 
for a discussion of the relative merits of each mechanism. In 
this study we are concerned with the effect of a PMZ and the 
13 C(a, n) 16 on light-element nucleosynthesis. The neutrons 
from this reaction are important for the prod uction of e.g., 
19 F l|Forestini et al.lll992l ; ILugaro et al.ll2004l ) The method 



i „ ,Y ' ' ' 9 r, , . ' ,' on,, , ,% . ,„ m , — : — we use to include a PMZ has been described in Lugaro et al.l 

St " MS(P ' 7 > ^ M ^ ' 7) ^ 7J^ ele » and is similar to that used by Ic^vftMmlavi 



26 Al 9,m are the ground and meta-stable state of 26 Al re- 
spectively), the 26 Mg(p,7) 27 Al and the 27 Al( p, 7) 28 Si and 
27 Al(p , a) 24 Mg reaction rates are taken from llliadis et all 
(|200lT ) . Other important updates include the most re- 
cent evaluation of t he 14 N(p,7) 15 CNO rate from 
2006), and the revised triple-a rate from 



Bemmerer et al.l ._ L 
Fvnbo et al l (|2005h . Many of the rate s are now from 



the NACRE compilation (lAngulo et all [l999) including 
the following CNO reactions: 12 C(p, 7) 13 N, 13 C(p, 7 ) 14 N, 
16 0(p, 7 ) 17 F, 18 0(p, 7 ) 19 F and ls O(p,«) 15 N rates. The 
O fa 7) 18 F and 17 Q(p , a) 14 N rates are from the more re- 
cent [ChafoeiaD l|20071 ). 

The neutron source reactions have also been updated in- 
cluding the 13 C(a,n) 16 rate wh ich is taken fr o m NA CRE, 



whereas previously we used the Denker et all |l995) rate. 
Previously we used the iKaeppeler et al.1 (|l994h rates for the 
22 Ne(a, n) 25 M g and 22 Ne(a, 7) 26 M g reactions, now we use 
the rates from IKarakas et al.l (|2006l ) . 

Some of the important /3-decay rates have also been 
update d or corrected. For example, the REACLIB database 
used in IKarakas fc Lattanziol (|2007f) had an incorrect decay 
rate for 60 Fe of 0.22 Myr; this has been corrected to 1.5 Myr 
in the current JINA REACLIB library. 



2.2 The inclusion of a partial mixing zone 

In AGB stars there are two important neutron produc- 
ing reactions. The 22 Ne(a, n) 25 Mg reaction operates dur- 
ing the convective thermal pulses when T > 300 x 10 6 K. 
This has been suggested to be the dominant neutron source 
in intermediate-mass AGB stars, whereas these tempera- 
tures are reached only in the last few TPs of lower mass 
(M < 3Mq) stars. The other potential source of neutrons 
in AGB stars is the 13 C(a, n) ls O reaction, which operates at 
lower temperatures (T > 90 x 10 6 ) than the 22 Ne source. Ob- 
servational and theoretical evidence suggests thi s is the dom- 
inant neutron source in low-mass AGB stars l|Smith et al.l 
ll987l ; lGallino et al.lll998h . 

To operate efficiently the 13 C(a, n) le O reaction requires 



(|2000|). We include a PMZ of constant mass of O.OO2M at 
the deepest extent of each third dredge-up episode for the 
3M , Z = 0.02, 0.008, and 0.004 models, and for the 2M , 
Z = 0.0001 model. 

Here we only examine the effect of the PMZ on low-mass 
AGB models. This is because we add the partially-mixed 
zone in the post-processing step, so any potential feedback 
on the structure of the star is ignored. In intermediate-mass 
models of low-metallicity (M > 3Mq at Z < 10~ 4 ) the 
temperature at the base of the convective envelope dur- 
ing dredge-up may become hot enou gh for proton-captures 
l|Gorielv fc Siess 2004; Hcrwig 2004a|). The energy produced 
by these hot dredge-ups may effect the st ructure of the star, 
by increasing the depth of dredge -up dHerwia 2004at). o r 
by terminating the AGB altogether (|Woodward et al.ll2008l ) . 
This situation could arise if the ingestion of protons leads 
to an hydrog en flame that produces enough energy to eject 
the envelope (I Woodward et ai1l2008l '). Consequences of pro- 
ton ingestion on the nucleosynthesis is largely unknown but 
could include the inhi bition of formation of the 13 C pocket 
l|Gorielv fc Siessll2004l ). 

The models we consider in this study do not fall into the 
region of the mass-Z-plane wh ere hot dredge-ups ar e pre- 
dicted to occur (see Fig. 4 from lGorielv fc Si ess 2004). The 
3M , Z — 0.004 model is the only candidate for hot dredge- 
ups and in this model the peak temperature is < 30 x 10 6 K 
during a few dredge-up episodes, but accompanied by low 
densities (p < lg/cm 3 ). It is unlikely that proton captures 
occur during the TDU under these conditions, although fur- 
ther study into the inclusion of a PMZ on the structure of 
the 3Mq, Z = 0.004 model is required. Note that we do not 
include a PMZ into the 3M , Z = 0.0001 model, as this 
would suffer hot dredge-ups, where the peak temperature at 
the base of the envelope during dredge-up is ~ 60 x 10 6 K. 

Even if proton ingestion is unlikely to occur in our low- 
mass AGB models, these studies do point out the inconsis- 
tency of adding a PMZ into the post-processing calculation. 
For this reason the results presented in Section 14.21 should 
be treated with some caution. The trends described in this 
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Table 1. Grids of stellar masses for each Z, noting if the models 
experience the core He-flash (CHe), the third dredge-up (TDU), 
and hot bottom burning (HBB). The second line lists the mass- 
loss law used on the AGB, including the parameter r) if the 
Reimers mass-loss formula was used. 



1V± dot) 


Z = 0.02 


Zj — U.UUO 


— u.uu^ 


Z = 10 — 4 


1.0 


CHe 


CHe 


CHe 


CHe, TDU 




VW93 


VW93 


VW93 


VW93 


1.25 


CHe 


CHe 


CHe 


CHe, TDU 




VW93 


VW93 


VW93 


VW93 


1.5 


CHe 


CHe 


CHe, TDU 


CHe, TDU 




VW93 


VW93 


VW93 


VW93 


1.75 


CHe 


CHe, TDU 


CHe, TDU 


CHe, TDU 




VW93 


VW93 


VW93 


VW93 


1.9 


CHe 


CHe, TDU 


CHe, TDU 


CHe, TDU 




VW93 


VW93 


VW93 


VW93 


2.0(2.1) 


CHe 


CHe, TDU 


CHe, TDU 


TDU 




VW93 


VW93 


VW93 


VW93 


2.25 


CHe,TDU 


TDU 


TDU 


TDU 




VW93 


VW93 


VW93 


VW93 


2.5 


TDU 


TDU 


TDU 


TDU 




VW93 


VW93 


VW93 


VW93 


3.0 


TDU 


TDU 


TDU 


TDU, HBB 




VW93 


VW93 


VW93 


R75,ry = 5 


3.5 


TDU 


TDU 


TDU 


TDU, HBB 




VW93 


VW93 


VW93 


R75,?? = 7 


4.0 


TDU 


TDU, HBB 


TDU, HBB 


TDU, HBB 




VW93 


VW93 


VW93 


R75,7? = 7 


4.5 


TDU, HBB 


TDU, HBB 


TDU, HBB 


TDU, HBB 




VW93 


VW93 


VW93 


R75,?7 = 10 


5.0 


TDU, HBB 


TDUHBB 


TDU, HBB 


TDU, HBB 




VW93 


VW93 


VW93 


R75,7? = 10 


5.5 


TDU, HBB 


TDUHBB 


TDU, HBB 


TDU, HBB 




VW93 


VW93 


VW93 


R75,?7 = 10 


6.0 


TDU, HBB 


TDU, HBB 


TDU, HBB 


TDU, HBB 




VW93 


VW93 


VW93 


R75,7? = 10 


6.5 


TDU, HBB 
VW93 









paper are likely to be qualitatively accurate but the actual 
extent of the effect of a PMZ may be different to that de- 
scribed here. Further detailed studies into the effects of a 
partial mixing zone on AGB nucleosynthesis is required, and 
will be the subject of a future investigation. 



3 THE AGB MODELS 

The grids of AGB models used to compute the updated 
yields are provided in Table [1] In this table it is noted if the 
model experienced the core He-flash, TDU and/or HBB, and 
the mass loss formulation employed on the AGB. Most of the 
models list ed in Table [J have been previous l y described in 
detai l (e.g.. lKarakas fc Lattan zio 2003, 2007J; Karakas et al.l 
2009). Here we restrict the discussion to the new AGB mod- 
els, which are the 4.5 and 5.5M , Z = 0.02, 0.008 and 0.004, 
the 2.1M , Z = 0.004, and the 1.0, 1.5, 1.9Mq, Z = 0.0001 
models with VW93 mass loss. The new intermediate-mass 
Z = 0.0001 AGB models of 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, and 
6Mq are computed using R75 as described in Section [2] 



The parameter ?7r used for each model is listed in Table [T] 
These parameters were chosen as follows. Little is known 
about mass loss from low-metallicity AGB stars of interme- 
diate mass. However, detailed evolution and nucleosynthesis 
models of such stars indicate that the total C-r-N+O content 
of the envelope quickly reaches similar levels found in AGB 
stars of much higher metallicity o wing to a very efficient 
TDU (see in particular, §2.1 from iHerwial l2004bh . Hence 
there is no good reason to suspect that the mass loss from 
these stars should be significantly lower than, say, AGB stars 
of Z — 0.004. For this reason, we choose the »7r parameters 
such that the final number of thermal pulses is similar to 
that found in Z — 0.004 models of the same mass. We note 
that these choices are arbitrary, and changes in ??r will have 
a strong impact on the chemical yields. 

In Table [2] we present some structural details from the 
new models including the initial mass and metallicity, the 
final core and envelope mass (M COIC and M auv , respectively), 
the number of TPs computed, and the total amount of mat- 
ter dredged up into the envelope during the TP- AGB phase 
(Mj°* d ). All masses are in solar units. There are two entries 
for the 4Mq, Z = 0.0001 case. The first lin e shows the re- 
sults c omputed with the code that includes iFerguson et al.l 
|2005l ) low-temperature opacities, the second with the older 
version used to compute the models listed above that in 
Table [5] The most important parameters, from a nucle- 
osynthetic point of view, are the total number of thermal 
pulses as this determines how much material from the He- 
shell is mixed into the envelope, and the efficiency of HBB. 
These numbers are similar, although the model computed 
with the older code ended with a higher final envelope mass 
and could, in principle, experience at least one more TP. 
The e ffect of the (scaled solar composition) Fer guson et al.l 
l|2005l ) low-temperature opacities on the evolution is small at 
such low metallicities, and is the main reason for the consis- 
tency between the two 4Mq models. The use of carbon-rich 
low temperature opacities may have a much larger impact, 
owing to the fact that all of the Z — 0.0001 models, even 
the IMq and 6Mq models, become carbon rich with a final 
C/O > 1. 

We provide an electronic on-line table with details of 
each new stellar model as a function of TP number; (Ta- 
ble Al). This data a r e sim ilar to the tables presented in 
iKarakas fc Lattanzid l|2007h but with the inclusion of the 
radius, bolometric luminosity, and effective temperature. 
These tables include, for each TP, the core mass (Mcore), 
the maximum extent of the flash-driven convective region 
(M cs h), the duration of flash-driven convection (t cs h), the 
amount of matter dredged into the envelope after that 
pulse (Ddredge), the thir d dredge-up efficien c y par ameter 
(A), and Ldup defined bv lGorielv fc Mowfavil i|200Ch to be 
AMdredge/Afcsh- The TDU efficiency parameter, A, is usu- 
ally defined according to A = AMdredge/AMh, where AMh 
is the amount by which the core mass has grown between 
the present and previous TPs. 

Further, we provide the maximum temperature in the 
He-shell (THeshell), the maximum temperature at the base 
of the convective envelope during the previous interpulse 
period (Tbce), the maximum temperature in the H-shell 
(THshell), the interpulse period, the total mass at the be- 
ginning of the TP, the maximum radiated luminosity during 
the previous interpulse period (MaxL), the maximum lumi- 
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Table 2. Details of the new AGB models. 



1V-L CLOO 


z 


M 


M 


No. TP 


M tot 

drcd 


4.5 


0.02 


0.853 


0.692 


21 


0.053 


5.5 


0.02 


0.900 


0.800 


36 


0.075 


1.5 


0.008 


0.861 


0.670 


38 


0.141 


5.5 


0.008 


0.907 


0.864 


56 


0.143 


2.1 


0.004 


0.650 


0.090 


20 


0.058 


1.5 


0.004 


0.873 


1.160 


50 


0.177 


5.5 


0.004 


0.932 


0.484 


71 


0.154 


1.0 


0.0001 


0.726 


8.5(-3) 


26 


2.6(-3) 


1.5 


0.0001 


0.662 


0.022 


18 


0.06 


1.9 


0.0001 


0.682 


0.029 


21 


0.170 


3.0 


0.0001 


0.812 


0.120 


20 


0.106 


3.5 


0.0001 


0.854 


0.136 


27 


0.103 


4.0 


0.0001 


0.872 


0.133 


37 


0.128 






0.876 


0.477 


36 


0.115 


1.5 


0.0001 


0.898 


0.085 


11 


0.109 


5.0 


0.0001 


0.929 


0.123 


56 


0.122 


5.5 


0.0001 


0.966 


0.085 


77 


0.124 


6.0 


0.0001 


1.008 


0.183 


109 


0.127 



nosity generated by the TP (MaxLHe) , the maximum radius 
during the previous interpulse period (maxR) , the bolomet- 
ric luminosity (Afbol) , and finally, the effective temperature 
(Teff). All units are in solar units with the exception of 
temperatures, that are in kelvin, and all times which are in 
years. 

We briefly comment on the 1M , Z = 0.0001 model. 
This model experienced 26 TPs and had a final core mass 
of 0.72Mq, higher than found in metal- rich models of the 
same initial mass (e.g., O.56M at Z = 0.02). The IMq, 
Z = 0.0001 experienced shallow TDU after the 5 th and 6 th 
TPs (the core mass was » 0.56M Q at the 6 th TP). The small 
amount of TDU was enough to make the star carbon rich, 
where C/O > 1. The addition of primary 12 C into the enve- 
lope had a dramatic effect on the evolution of the star. This 
is because the extra 12 C changed the rate of nuclear burn- 
ing in the H-shell. The additional energy caused a significant 
shortening of the interpulse period from ~440,000 years at 
the 6 th TP to 240,000 years, followed by a steady decrease 
to 45,000 years, and there was no further TDU. S imilar be- 
haviour was reported by IStancliffe fc Glebbeekl (|200&f ) for 
the same mass and metallicity. 

For the remainder of this section we compare the R75 
model s to the VW93 models from iKarakas fc Lattanziol 
(|2007h . We will use the 3 and 6M , Z = 0.0001 mod- 
els as examples. The 3M model with VW93 had 40 TP 
and dredged a total of 0.242Mq into the envelope. This 
model also experienced HBB with a peak temperature of 
71 x 10 6 K. The final core mass was 0.82Mq. In compari- 
son, the model with R75 mass loss had 20 TPs and dredged 
0.106Mq into the envelope. This model only had mild HBB, 
with a peak temperature of ~ 40 x 10 6 K. The final core 
mass was 0.8 IMq , probably the only similarity to the VW93 
model. 

In contrast, both 6Mq, Z — 0.0001 models had about 
the same number of TPs: 106 in the VW93 case compared to 
109 in the case of the R75 model. The peak temperature at 
the base of the envelope was also similar at 104 x 10 6 K, along 




6.23xl0 7 6.24xl0 7 6.25xl0 7 

Time (in years) 



Figure 1. The temporal evolution during the TP- AGB of the 
temperature at the base of the convective envelope for 6Mq Z = 
0.0001 models with R75 (top panel) and VW93 (lower panel) mass 
loss. 

with the total amount of matter dredged into the envelope, 
(0.114M Q for the VW93 model compared to 0.127M Q ). The 
real difference between the models is that the R75 model was 
evolved to a small envelope mass of only 0.13Mq, which led 
to the cessation of HBB. This can be seen in Figure [1] where 
we illustrate the evolution of the temperature at the base of 
the convective envelope during the TP- AGB for the R75 (top 
panel) and VW93 (lower panel) models. In the top panel the 
temperature drops dramatically to below 10 x 10 6 K during 
the last few TPs. This is caused by the rapid reduction in 
the envelope mass. In comparison the VW93 model (top 
panel) only loses ~ 0.04Mq during the AGB and we expect 
that this model may experience at least another ~ 100 TPs , 
setting the yields presented in lKarakas fc Lattanziol (|2007l ) 
as lower limits. 



4 STELLAR YIELDS 

To compute the yields we integrate the mass lost from the 
model star during the entire stellar lifetime according to 

M i = f [X(i)-X (k)]^-dt, (3) 

where Mi is the yield of species i (in solar masses), dM/dt 
is the current mass-loss rate, X(i) and Xo(i) refer to the 
current and initial mass fraction of species i, and r is the 
total lifetime of the stellar model. The yield can be negative, 
in the case where the element is destroyed, and positive if 
it is produced. We also present the total amount of i (in 
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Mq) expelled into the ISM, noting that this value is always 
positive. 

The stellar yields for the models listed in Table [T] are 
available as electronic on-line tables. We provide one table 
for each metallicity (Tables A2 to A5). in the Appendix for 
an example. Each table contains the following: 1) the nu- 
clear species, i, 2) the atomic mass, A(i), 3) the net stellar 
yield defined above, 4) the amount of species i in the wind 
lost from the star, mass(i)i ost , 5) the amount of i that would 
have initially been present in the wind, mass(i)o. The quan- 
tity mass(i)o is the mass expelled during the stellar lifetime 
multiplied by the initial mass fraction. We next include 6) 
the average mass fraction of i in the wind, (X(i)), 7) the 
initial mass fraction Xo(i), and 8) the production factor / 
defined by log 10 [(X (i)) / Xo(i)]. All yields are in solar masses 
and abundances in mass f raction. These yields a r e pres ented 
in the same format as in iKarakas fc Lattanziol (I2007T ) . We 
note that there are two entries for 26 Al: the first shows the 
yields for the longer- lived ground state 26 A1 S which has a 
half life of 7.17 x 10 5 years and is denoted by the symbol 
"al-6" in the on-line tables; "al*6" refers to the short-lived 
meta-stable state. The yields of "al-6" should be added to 
that of 26 Mg for chemical evolution studies. 



In lKarakas fc Lattanziol l|2007i ) we discussed the inclu- 
sion of of synthetic TPs to account for TPs not computed in 
detail. This was done because some of the models had rea- 
sonably large envelope masses at the end of the computation 
and could, in principle, experience further TPs and dredge- 
up. We do not repeat that discussion here but we also do not 
include the contribution of these synthetic TPs for the fol- 
lowing reason. We would need to make an assumption about 
the efficiency of the third dredge-up parameter at small en- 
velope mass, and this is a unknown. There is evidence to 
suggest that A decreases with decreasing envelope m asses 
for low-mass l|Straniero et al . 1997; Karakas et alj2002h and 
intermediate-mass AGB stars l|Vassiliadis fc Wood! fl993l ). 
The new intermediate-mass Z — 0.0001 AGB models pre- 
sented here were evolved to small envelope masses (~ O.IMq 
in some cases). These models do not in general show any 
decrease in A with decreasing envelope mass. However, it 
should be pointed that low-metallicity models tend to ex- 
perience mor e efficient mixing than metal-rich AGB stellar 
mode ls (e.g., iBoothrovd fc Sackmannl Fl988l ; IKarakas et alj 
l2002h . 

O ne last difference between t he yields presented here 
and in IKarakas fc Lattanz io (2007j) is that we do not specif- 
ically provide yields for planetary nebulae (PNe). That is 
because column 6 of the on-line yield tables provides the 
average mass fraction of the wind. This value is weighted 
toward the composition of the envelope during the last few 
TPs, because that is when most of the mass is lost from 
the star. Hence this value is suitable for comparison to PNe 
abundances. The models computed using the R75 mass-loss 
formula have not had their mass loss weighted toward the 
final few TPs, as they would have if we were to use the 
VW93 prescription. However few, if any, PNe should exist 
that have evolved from intermediate-mass Z — 0.0001 AGB 
stars. The few PNe found in the Galactic Halo likely evolved 
from stars with initial masses closer to ~ IMq. 



4.1 Comparison to previous work 

In Figures [2] to [5] we show the weight ed yields of 12 C, 14 N, 
17 Q, a nd 19 F from this study and from IKarakas fc Lattanziol 
(|2007h for the Z = 0.02, 0.008, 0, 004, and Z = 0.0001 mod- 
els, respectively. In each plot, the yields have been weighted 
by th e initial mass function (IMF) from iKroupa et al.l 
J1993J). Note for this comparison that we are using the yields 



Karakas fc Lattanziol l|2007l ) with no contribution from 



from 

synthetic TPs. These isotopes were chosen to be representa- 
tive of light element nucleosynthesis in AGB stars of various 



At Z 



= 0.02, the new yields of 12 C, 17 0, and 
19 F are reasonably cons istent with the yields from 
IKarakas fc Lattanziol (|2007h . We see small reductions in the 
yield of 12 C for M > 2.5Mq, and small increases in the 
yields of 19 F for M > 4M . The yields of 14 N are system- 
atically lower for masses less than 4Mq , whereas above this 
we can see the effect of a finer mass grid on the shape of 
the yield as a function of stellar mass. For the Z — 0.008 
and 0.004 models (Figures [3] and [4| we see a similar trend as 
seen for Z = 0.02, although the new yields of 17 and 19 F 
are generally larger. The 12 C and 14 N yields are remarkably 
consistent, indicating that the change of reaction rates or 
initial abundances has not had a significant impact on the 
yields of these important isotopes. 

The new Z = 0.0001 model yields are consistent with 
the old yields for 12 C and for 19 F. The main difference for 
fluorine is that the peak production is shifted to 2.5Mq from 
2.25Mq, as a conseq uence of using the lo wer 19 F(a,p) 22 Ne 
destruction rate from lUgalde et all (|2008l ) . The main differ- 
ence between the yields in Figure [5] are for 14 N and 17 0. In 
both cases, the new yields produce less of each isotope at 
a given mass for M > 2.5Mq. The differences observed at 
2.5Mq are surprising, since we are using the same structure 
model. The changes in the yields at this mass are a conse- 
quence of the new reaction rates (and in particular the new 
14 N(p, 7) 15 and 17 +p rates). However the difference for 
14 N is most apparent at 3Mq , where the variation is a factor 
of ~ 650. If we assume that the rates have introduced at least 
a factor of two variation, then the change of structure model 
(using the R75 model instead of the VW93) is still by far the 
most significant factor. For 17 0, the new intermediate-mass 
Z = 0.0001 models result in smaller yields for M > 3Mq, 
with the biggest variation at 4 — 5Mq . 

In Figure[S]we show the weighted yields of 23 Na for the 
Z = 0.02, 0.008, 0.004, and 0.0001 models, again compar- 
ing th e new yields to the yields from IKarakas fc Lattanziol 
(2007). It is for this element that we see the largest differ- 
ences as a result of using the new reaction rates. The new 
yields are systematically lower for all stellar models com- 
puted. Table [3] shows that the variation increases with in- 
creasing mass, and at a given mass, with decreasing metal- 
licity. That is, models with HBB show the largest impact 
of using the n ew 23 Na(p, 7) 24 Mg and 23 Na(p, a) 20 Ne reac- 



tion rates from Hale et al. (2004), which are faster than the 



lEl Eid fc Champagne! fll995) rates that we used previously. 
This has resulted in higher rates of 23 Na destruction and 
dramatically lower yields of sod ium for M > 3Mq for all 
metallicities. Ilzzard et all l|2007h found a similar result us- 
ing both detailed an d synthetic AGB models an d a similar 
set of reaction rates. IVentura fc D'Antonal |2005l ) discussed 
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Figure 2. Weighted yields of 12 C, 14 N, 17 0, and 19 F as a function of the initial mass for the Z = 0.02 models. The solid line and open 
circles show results for the updated yields; the dashed line and filled circles show results from Karakas & Lattanzio (2007). 
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Figure 3. Weighted yields of 12 C, 14 N, 17 0, and 19 F as a function of the initial mass for the Z = 0.008 models. Symbols are the same 
as in Figure [2] 
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the variation of 23 Na production as a re sult of uncertainties 
in rea ction rates. They compared the ICaughlan fc Fowlerl 
(1988) rates to the NACRE database, and found similarly 
large differences as reported here. 



4.2 The effect of a partial mixing zone 

In this section we consider the effect of a partially mixed 
zone inserted at the deepest extent of each third dredge-up 
episode on the stellar yields. In Table [4] we show the yields 
for selected isotopes from the 3Mq models. The yields for 
all species for each model with a PMZ is included as an 
on-line data table (Table A6). The effect of a PMZ on the 
nucleosynthesis of a 2Mq, model of Z = 0.0 001 is discussed 
in Sec tion 14.31 in comparison to results from ICristallo et al.l 
(2009 ). In each model, a PMZ of 0.002M Q has been inserted. 
Note that in the computations the same initial abundances 
and reaction rate network are used. The only variation is 
the inclusion of a partially mixed zone at the deepest extent 
of each dredge-up episode. 

The extra 14 N in the intershell results in higher yields 
of 22 Ne, and hence also 23 Na. The production of 23 Na oc- 
curs via the 22 Ne(n, 7) 23 Ne(/3~ 1 ) 23 Na seq uence of reactions 
ijGorielv fc Mowlavj Eoool : Irlerwid [2004bh . and is the main 
explanation for the increase in this element in Table [4] The 
neutrons from the 13 C(a,n) 16 reaction are also responsi- 
ble for enrichments in other isotopes, including 19 F, 21 Ne, 
25 Mg, 26 Mg, 26 Al, 27 A1, 29 Si, 30 Si, 31 P, and 60 Fe. 

The isotopes of 31 P and 60 Fe are only produced by 
neutron capture in AGB stars. The production of the ra- 
dioactive 60 Fe requires high neutron densities to overcome 
the branching at 59 Fe. Hence 60 Fe production is primar- 
ily via the 22 Ne(a, n) 25 Mg source during TPs. Hence the 
3Mq, Z = 0.02 model without a PMZ produces very lit- 
tle P or 60 Fe, as it is not hot enough to sustain the 
22 Ne(a,n) 25 Mg source during TPs, and without a PMZ, 
does not have enough 13 C present to produce any neutron 
capture elements. The increase of 60 Fe caused by the PMZ 
is the result of the extra 22 Ne in the intershell (from 14 N) 
causing a slightly higher activation of the 22 Ne(a, n) 25 Mg 
reaction. The difference in the production of 31 P and 60 Fe 
is lessened in the lower metallicity models, mostly as a con- 
sequence of partial activation of the 22 Ne neutron source in 
models without a PMZ. 



4.3 Comparison to other authors 

In this section we compare the structure and yields from 
the 2 M0, Z = 0.0 01 m odel to the comparative model 
from ICristallo et al.l (|2009h . For the sake of the compari- 
son we compute a 2Mq, Z = 0.0001 model with a PMZ 
of 0.002M Q , which results in a 13 C pocket that is ~ 10% 
of the mass of the in t ershe ll. This is because all the mod- 
els in ICristallo et afl ()2009h have 13 C pockets that formed 
as a consequence of the algorithm used to find the convec- 
tive border following a TP. We assumed an exponentially- 
decaying proton profile in our PMZ models and this results 
in a 13 C pocket that is s imilar in shape to the [3 = 0.1 case 
in ICristallo et al.1 (|2009l . see the top panel of Fig. 4). The 
P — 0.1 case is their adopted value in the full detailed nu- 
cleosynthesis computations. An important difference is the 



mass extent of the p ocket. From the top panel of Fig. 4 in 
ICristallo et al.l (|2009l ). we estimat^l that the extent of the 
13 C pocket (where the 13 C abundance is larger than the 14 N 
abundance) to be » 0.0015Mq, whereas in our model we 
obtain a value about 1.5 times smaller, at ~ 9 x 10~ 4 Mq. 
This mass difference has important implications for the nu- 
cleosynthesis as discussed below. 

The main difference between the models is the 
treatment of convective borders, and the use of car- 
bon and nitrogen-ri ch low-temperature opacities in the 
ICristallo et al.l |2009Ti model. The effec t of C-rich opa cities 
on AGB evolution was first discussed bv lMarigol £2002), who 
noted that once C/O > l the star cools and mass loss in- 
creases, shortening the TP- AGB lifetime. This has the effect 
of reducing the total number of TPs (Cristallo et al. find 15, 
we have 26), and the total amount of mass dredged up into 
the envelope (our model dredges up 30% more He-shell ma- 
terial) . 

The structural differences between the ICristallo et al.l 

(2009) model and ours is reflected in the stellar yields pro- 
vided in Table [5] Cristallo et al. also used different initial 
abundances, s caling their Z = 0.0 001 model to the solar 
abundances of lAsplund et al.l j|2005l ) (Z Q = 0.0138, smaller 
than our Zq = 0.02). The differences caused by this are 
hard to asses, as they did not provide their initial abun- 
dances (ours are given in column 7 in the on-line yield 
data tables). Even so, we find relatively large variations 
for all isotopes, even those not directly produced by He- 
burning such as 14 N. In these low-metallicity models, most 
of the 14 N is produced during the AGB, as a consequence of 
the dredge-up of material processed by CNO cycling in the 
H-shell. Our yields of He-burning products including 12 C, 
ie O, 19 F, are larger as expected from our model dredging 
up more He-intershell material. The increase in the surface 
abundance of 16 O results in a final [O/Fe] ~ 1.0, and the 
planetary nebulae formed from such an object would show 
the dredge-up of oxygen. The metallicity dependent enrich- 
ment in O predicted by our models has been reported by 
iMagrini fc Goncalvesl (|2009h for planetary nebulae in the 
galaxy IC 10. 

Back to the comparison to Cri stallo et al. I (|2009l) . while 
our yields of He-burning products are larger, the products 
of neutron captures (e.g., 31 P) along with 26 Al and 27 A1 are 
smaller. This is a consequence of our smaller 13 C pocket (at 
least as measured at the 2 nd TDU). The variations range 
from only -14% for 26 Al up to ~ -90% for 60 Fe. These 
variations may also reflect different choices of reaction rates 
(e.g., for Al). One interesting difference is that our 2Mq 
model with a PMZ produces a positive yield of 15 N, whereas 
both the model without a PMZ and the ICristallo et al.l 
(2009) model have negative yields, indicating a net destruc- 
tion. The positive yield indicates that some 15 N is produced 
during TPs by the 18 0(p, a) ls N reacti on rate, which is par t 
of the fluorine production chain (see iLugaro et alj |2004| ) . 
The difference in the 15 N surface abundance between the 
models with and without a PMZ is ~ 50%, and most of this 
extra 15 N is produced during the first few TPs when the 



2 these mass estimates are smaller than the effective 13 C pocket 
size denned in Cristallo et al. (2009). In either case our 13 C pock- 
ets are somewhat smaller. 
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Table 3. Yields (in Mq) of 23 Na from models with M > 3Mq . The old yield refers to the yields from Karakas & Lattanzio ( 20 0?]). For the 
Z = ,0001 models, we compare the sodium yields from the R75 models presented here to the VW93 models from [Karakas &~Lattanziol 
(2007). For the 6M0, Z = 0.0001 there are two entries: the first shows the new and old yields from the VW93 model, and the second 
line the new yield for the R75 model. All yields are expressed in the form n(m) = n x 10 m . 



Initial mass 


Z 


New 23 


Na yield 


Old 23 Na yield 


Factor difference 


3.0 


0.02 


5.0942( 


-5) 


7.0317(-5) 


1.38 


3.5 


0.02 


6.5188( 


-5) 


9.3885(-5) 


1.44 


4.0 


0.02 


5.5071( 


-5) 


6.7809(-5) 


1.23 


5.0 


0.02 


8.9715( 


-5) 


1.3889(-4) 


1.55 


6.0 


0.02 


1.3131( 


-4) 


1.0715(-3) 


8.16 


6.5 


0.02 


1.4863( 


-4) 


1.0026(-3) 


6.75 


3.0 


0.008 


7.3897( 


-5) 


1.4416(-4) 


1.95 


3.5 


0.008 


3.2963( 


-5) 


6.2545(-5) 


1.90 


4.0 


0.008 


2.5255( 


-5) 


4.9604(-5) 


1.96 


5.0 


0.008 


6.4444( 


-5) 


1.7090(-3) 


26.5 


6.0 


0.008 


4.4473( 


-5) 


8.7833(-4) 


19.7 


3.0 


0.004 


4.3768( 


-5) 


7.2967(-5) 


1.67 


3.5 


0.004 


1.7585i 


-5) 


2.9209(-5) 


1.66 


1.0 


0.004 


2.5553( 


-5) 


1.5553(-4) 


6.09 


5.0 


0.004 


6.8179( 


-5) 


1.6316(-3) 


23.9 


6.0 


0.004 


1.0933( 


-5) 


4.4738(-4) 


40.9 


3.0 


0.0001 


1.3794( 


-5) 


1.7719(-3) 


128 


3.5 


0.0001 


1.2139( 


-5) 


2.6254(-3) 


216 


1.0 


0.0001 


1.7816( 


-5) 


2.5865(-3) 


145 


5.0 


0.0001 


1.0310( 


-5) 


1.8136(-3) 


176 


6.0 


0.0001 


7.0570( 


-6) 


7.6308(-5) 


10.8 


6.0 


0.0001 


1.0482( 


-5) 
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Table 4. Yields for selected isotopes from the 3Mq models. The first line shows the yield from the model with no PMZ, the second line 
the yield from the PMZ model, and the third line the percentage difference between the yields. 



z 


19 p 




22 Ne 




23 Na 




25 Mg 




26 Mg 




30 Si 




31p 




60 Fe 




0.02 


3.4087(- 


-6) 


1.8092(- 


-3) 


5.0942(- 


-5) 


1.0767(- 


-5) 


1.3568(- 


-5) 


7.7026(- 


-7) 


8.0023(- 


-8) 


9.5257(- 


-9) 




4.1292(- 


-6) 


3.4512(- 


-3) 


9.2811(- 


-5) 


2.1733(- 


-5) 


2.5835(- 


-5) 


3.1746(- 


-6) 


1.6082(- 


-6) 


2.2711(- 


-8) 




21.1% 




90.8% 




82.2% 




102% 




90.4% 




312% 




1910% 




138% 




0.008 


1.1481(- 


-5) 


4.2290(- 


-3) 


7.3897(- 


-5) 


6.7804(- 


-5) 


7.3233(- 


-5) 


2.5349(- 


-6) 


2.7079(- 


-7) 


3.4270(- 


-7) 




1.3700(- 


-5) 


7.5472(- 


-3) 


1.6357(- 


-4) 


1.0954(- 


-4) 


1.3010(- 


-4) 


3.1024(- 


-6) 


2.1346(- 


-6) 


4.0355(- 


-7) 




19.3% 




78.5% 




121% 




61.6% 




77.6% 




22.4% 




688% 




17.8% 




0.004 


1.1311(- 


-5) 


3.4789(- 


-3) 


4.3768(- 


-5) 


1.0715(- 


-4) 


1.4453(- 


-4) 


2.9138(- 


-6) 


4.1899(- 


-7) 


1.4891(- 


-6) 




1.4083(- 


-5) 


6.9083(- 


-3) 


1.0657(- 


-4) 


1.8292(- 


-4) 


2.7583(- 


-4) 


3.8189(- 


-6) 


1.9773(- 


-6) 


1.7557(- 


-6) 




24.5% 




98.6% 




143% 




70.7% 




90.8% 




31.1% 




372% 




18.0% 





Table 5. Yields for selected isotopes from the 2Mq Z = 0.0001 models. The first line shows our model without a PMZ, the second line 
our model with a PMZ of 0.002Mq, and the third line the yields from the Cristallo et al. model (C09). We also show the percentage 
difference between line two (PMZ model) and the line three (C09 model). 



Model 


12 C 




14 N 




15 N 




ICQ 




19 p 




23 Na 




25 Mg 




30 Si 


31p 


60p e 




No PMZ 


3.57(- 


2) 


7.08(- 


5) 


-3.35(- 


-9) 


5.15(- 


-4) 


1.10(- 


-5) 


1.04(- 


-4) 


4.87(- 


-5) 


6.53(-8) 


1.39(-8) 


1.18( 


-8) 


PMZ 


3.27(- 


2) 


6.73(- 


5) 


1.07(— f 


*) 


9.52(- 


-4) 


1.30(- 


-5) 


2.26(- 


-4) 


7.64(- 


-5) 


9.14(-8) 


1.30(-7) 


7.04( 


-9) 


C09 


1.71(- 


2) 


3.40(- 


5) 


-2.07( 


-8) 


4.03(- 


-4) 


2.44(- 


-6) 


1.38(- 


-5) 


2.55(- 


-5) 


9.89(-8) 


3.53(-7) 


5.76( 


-8) 




91.2% 




98.0% 




152% 




136% 




433% 




154% 




200% 




-7.60% 


-63.2% 


-87.£ 


% 



temperature is not yet hot enough for efficient 15 N destruc- 
tion (or 19 F production). By the last TPs, the 15 N produced 
in the convective TP is rapidly destroyed in models with and 
without a PMZ. 

IStancliffe fc Jeffervl (120071 ) made a careful compar- 



ison between the 1.5M 







0.008 yields from 



stood by considering that the 22 Ne(a, n) 25 Mg reaction is 
margina lly activated during TPs at th is mass and metallic- 
ity, and IStancliffe &: Glebbeekl (|2008l ) used the faster rate 
from Kaeppeler et alj (|1994l ). That the yields for 23 Na are 
the same is a coincidence. Our model produces more 22 Ne 
and would produce more 23 Na, but we use the new faster 



Karakas fc Lattanziol ifeOOTfl t o their VW93 model, whereas Na P roton destruction rates 



Stancliffe fc Glebbeekl 1 20081 ) c ompared models of Z 



0.0001. The main difference that IStancliffe fc Jeffervl (|2007l ) 
found for the 1.5Mq, Z = 0.008 model was that their cal- 
culation experienced more efficient dredge-up, and hence 
positive net yields of 12 C, 22 Ne, and 25 Mg. In compari- 
son, our model had little dredge-up with only a total of 
1.47 x 10~ 3 Mf7i of He-shell ma terial mixed into the enve- 
lope. [Stanchffel&IIeff cry (200^) noted that the yields of iso- 
topes dependent on the efficiency of the first dredge-up (e.g., 
14 N) were reasonably consistent, given that their models had 
scaled -solar initial abundances and the lKarakas fc Lattanziol 
|2007l ) model did not. From our updated model with scaled- 
solar initial abundances, we find an agreement of ~ 5% be- 
tween the 14 N yields. 

IStancliffe fc Glebbeekl | 2008) have made a detailed 

comparison to the Z = 0.0001 models available in 

iKarakas fc Lattanziol (|2007l ); here we limit our discussion 
to the new 1.5Mq, Z — 0.0001 model. Our model has 
18 TPs compared to their 10, and dredged up ~ 52% more 
material from the He-intershell (R. Stancliffe, private com- 
munication). Hence our yields of 12 C, 14 N, and 22 Ne are 
larger by 9%, 64%, and 62% respectively. The Stancliffe 
& Glebbeck yields for 23 Na are almost identical, differ- 
ing by less than 1%, whereas their model produces almost 
100% more 25 Mg and 31 P. The last point can be under- 



Lastly, we compare the results from the new 
intermediate- mass Z = 0.0001 models to the models 
presented by IVentura" fc D'Antonal (I2009TI. These authors 
also p rovide a detailed comparison to lKarakas fc Lattanziol 
(2007) so here we comment on some of the major changes 
introduced by using the new R75 structure models and the 
large reduction in the Na yields. The new 4Mq, Z = 0.0001 
model has a similar number of TPs (37) to the 4Mq model 
computed bv lVentura fc D'Antonal |2009l ) (32), whereas our 
6M had 109 TPs compared to their 77 TPs. Further- 
more, owing to the differen t convective model used in the 
I Ventura fc D'Antonal |2009l ) calculations (the Full Spectrum 
of Turbulence compared to our use of the mixing-length the- 
ory with a = 1.75), the temperatures at the base of the 
envelope in our models were significantly lower. Our 4Mq 
model peaked at 88xl0 6 K compared to 92xl0 6 K; likewise 
our 6M n peaked at 104xl0 6 K compared to 127xl0 6 K in 
the I Ventura fc D'Antonal (|2009l ) model. Finally, one other 
major difference is evident in the efficiency of the third 
dredge-up where our models show A ~ 0.9 for M > 3Mq 
whereas in the Ventura & D'Antona models A varied from 
0.7 at 4M Q to 0.3 at 6M . 

These structural differences show up vividly in the 
chemical yields. All our Z = 0.0001 models have a fi- 
nal C/O > 1, whereas only the 6M , Z = 0.0001 
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I Ventura fc D'Antonal |2009l ) model has C/O > 1 and only as 
a consequence of O destruction. The final surface [O/Fe] of 
our 4Mq and 6Mq models are 1.19 and 1.14, respectively, 
compa red to 1.08 and —1.43 for the I Ventura fc D'Antonal 
(2009) models. The 4M yields for Na are the same, where 
the fin al surface [Na/Fe] abundance is 1.89 compared to 1.90 
for the I Ventura fc D'Antonal (2009) model. In contrast, the 
[Na/Fe] abundances from the 6Mq m odels are vastly differ- 
ent: o urs is 1.71 compared —0.20 from I Ventura fc D'Antonal 
(2009) model. This is almost entirely owing to their choice 
of convective model that results in much higher HBB tem- 
peratures and hence higher rates of 23 Na destruction. 



5 CONCLUSIONS 

In this paper we present new stellar yields of AGB models 
covering a range in initial mass from IMq to 6Mq , and ini- 
tial metallicity from Z = 0.02 to 0.0 001. These yields are 
an up date to the results presented in lKarakas fc Lattanziol 
(2007), using newer reaction rates and scaled-solar initial 
abundances for the Z = 0.008 and 0.004 models. We also 
present yields covering a finer mass grid than used previ- 
ously, and new intermediate-mass Z — 0.0001 AGB models 
between 3Mq and 6Mq. The new Z = 0.0001 models were 
evolved to low (~ O.IMq) envelope masses, with no reduc- 
tion in the third dredge-up parameter (A) observed. 

The main result of this paper is large reductions in the 
yields of 23 Na from models with HBB. This has implica- 
tions for the chemical evolution of galaxies and stellar sys- 
tems, and in particular for globular clusters. In contrast, 
the updated reaction network only results in small changes 
to the yields of 12 C, 14 N, and 19 F. The largest changes for 
these isotopes are caused by using new structure models at 
Z — 0.0001. These new structure models experience fewer 
TPs and TDU episodes, and this results in large reductions 
in the yields of 14 N and 17 0. 

We examine the effect of a partial mixing zone on the 
stellar yields of low-mass AGB models. Partial mixing zones 
are added into models of 2M Q , Z = 0.0001, and 3M , 
Z = 0.02, 0.008 and 0.004. The partial mixing zone results 
in a 13 C pocket in the top ~ l/10 th of the He-intershell 
and releases neutrons via the 13 C(a,n) le O reaction. The 
13 C pocket not only affects the production of 19 F but also 
other isotopes including 23 Na, the neutron-rich Mg and Si 
isotopes, 31 P, and 60 Fe. In particular, the yields of 31 P and 
60 Fe are the most affected by the introduction of the PMZ, 
where the yields of 31 P increase by factors of ~ 4 to 20 at 
3Mq depending on metallicity. Larger variations are found 
in the most metal-rich models. Because the PMZ is added 
into the post-processing step we ignore any feedback onto 
the stellar structure caused by the formation of the 13 C 
pocket. For this reason the yields should be treated with 
some caution. A full detailed study taking into account the 
feedback onto the stellar structure is necessary. 

The new result s are compared to the yields from 
iKarakas fc Lattanzid l|2007h . and to other recen t AGB 
yields in the lit e rature includi ng IStancliffe fc Jeffervl ( 
ICristallo et~all i|2009l) . and IVentura fc D'Antonal I 



2007), 



2009) 



The main result is that the structural differences between 
the calculations are still the dominant cause of variations 
between the nucleosynthesis predictions. In particular, the 



treatment of convection and the algorithm used to deter- 
mine convective borders are still important problems that 
need to be addressed in the future. 

One important piece of input physics missing from 
the models presented here is the inclusion of carbon and 
nitrogen-rich low-temperature opacities. While the approxi- 
mate treatment of molecular opacities we used in most mod- 
els does include some c ompositional depende nce, the low- 
temperature tables from iFerguson et al.l (120051 ) are for solar 
or scaled-solar mixtures only. The yields presented here and 
by many other authors show the complex nucleosynthesis 
that can occur during the TP-AGB (e.g., C/O > 1, C/N 
< 1). The composition of the envelope needs to be properly 
taken into account when determining the opacities of the 
outer layers. Such C and N-rich lo w-temperature opac i ty ta- 
bles are now becomin g available ( Lcdcrer & Armgerll2009l ; 
iHelling fc Lucasl 120091 ; iMarigo fc Aringerl 12009]) . and used 
in de tailed stellar-structure computations l|Cristallo et al.l 
2009). Stellar models using these opacities will need to be 
computed in the future for the size of the mass and metal- 
licity grid presented here. 
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APPENDIX A: EXAMPLES OF THE ON-LINE 
TABLES 

In Table [ATI we show the first few lines of the on-line table 
containing structural information from the new AGB mod- 
els. In Tables [H El IHIM] we show the first few lines of 
the on-line stellar yield tables for the Z = 0.02, Z = 0.008, 
Z = 0.004, and Z = 0.0001 models, respectively. An exam- 
ple of the on-line stellar yield for the models with partial 
mixing zones is shown in Table lAol 
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Table Al. The first few lines from Table Al. Each model entry starts with a header providing 
the initial mass and mctallicity of the model. To fit onto the page, we only show the first 8 
columns. 



#Minitial 


= 4.50 msun, 


Z = 0.0200 












#pulsc 




M csh 


^csh 


AM drcdgo 


A 


^dup 


Tkeshell 


1 

2 

3 : 


3.395468E-01 
3.407463E-01 
S.421246E-01 


2.890170E-03 
4.578114E-03 
4.609287E-03 


4.988529E+01 
3.937420E+01 
3.606038E+01 


0.000000E+00 
0.000000E+00 
1.928210E-04 


0.000000E+00 
0.000000E+00 
1.426807E-01 


0.000000E+00 
0.000000E+00 
4.183316E-02 


2.125615E+08 
2.376553E+08 
2.482608E+08 



Table A2. The first few lines from the Table A2. Each model entry begins with a header 
providing the initial mass and final mass (in solar units) along with the metallicity. 



# Minitial = 1.00 msun, Z = 0.0200, Mfinal = 0.564 msun 



Isotope i 


A 


yield 


mass(i) lost 


mass(i)o 


(X(i)) 


X0(i) 


/ 


9 a 


1 


8.1815390E-09 


3.2919932E-08 


2.4738393E-08 


7.5504431E-08 


5.6745975E-08 


1.2403737E-01 


n 


1 


-3.9796876E-43 


-3.9796876E-43 


0.0000000E+00 


-9.1224530E-43 


0.0000000E+00 


O.OOOOOOOE+00 


P 


1 


-8.1048310E-03 


2.9179156E-01 


2.9989639E-01 


6.6924673E-01 


6.8791503E-01 


-1.1948545E-02 



a g represents the sum of abundances from 64 Ni to Bi; an increase in g indicates that 
neutron-captures have occurred beyond the end of the network. 

Table A3. The first few lines from the Table A3. 



# Minitial = 1.00 msun, Z = 0.0080, Mfinal = 0.577 msun 



Isotope i 


A 


yield 


mass(i)i ost 


mass(i)o 


(X{i)) 


X0(i) 


/ 


g 


1 


3.3505394E-09 


1.2944828E-08 


9.5942889E-09 


3.0602433E-08 


2.2684230E-08 


1.3003190E-01 


n 


1 


0.0000000E+00 


O.OOOOOOOE+00 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


P 


1 


-8.7296963E-03 


3.0307159E-01 


3.1180128E-01 


7.1648133E-01 


7.3720652E-01 


-1.2384283E-02 



Table A4. The first few lines from the Table A4. 



# Minitial = 1.00 msun, Z = 0.0040, Mfinal = 0.610 msun 



Isotope i 


A 


yield 


mass(i)i ost 


mass(i)o 


(X(i)) 


X0(i) 


/ 


g 


1 


1.6834560E-09 


6.1053904E-09 


4.4219344E-09 


1.5654848E-08 


1.1339757E-08 


1.4004514E-01 


n 


1 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


O.OOOOOOOE+00 


P 


1 


-8.1304908E-03 


2.8376853E-01 


2.9189903E-01 


7.2761166E-01 


7.4855560E-01 


-1.2324413E-02 



Table A5. The first few lines from the Table A5. 



# Minitial = 1.00 msun, Z = 0.0001, Mfinal = 0.720 msun 



Isotope i 


A 


yield 


mass(i) lost 


mass(i)o 




X0(i) 


/ 


g 


1 


5.8190168E-09 


5.8983649E-09 


7.9347945E-11 


2.1065592E-08 


2.8343647E-10 


1.8711179E+00 


n 


1 


-3.2538150E-42 


-3.2538150E-42 


0.0000000E+00 


-1.1620968E-41 


0.0000000E+00 


0.0000000E+00 


P 


1 


-5.6426227E-03 


2.0989175E-01 


2.1553437E-01 


7.4961346E-01 


7.6990402E-01 


-1.1599216E-02 
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Table A6. The first few lines from the Table A6. 



# Minitial = 3.00 msun, Z = 0.0200, Mfinal = 0.682 msun, partial mixing zone = 2e-3 msun 



Isotope i 


A 


yield 


mass(i) lost 


mass(i)o 


(X(i)) 


X0(i) 


/ 


g 


1 


3.4544257E-06 


3.5859543E-06 


1.3152861E-07 


1.5470035E-06 


5.6745975E-08 


1.4355562E+00 


n 


1 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


0.0000000E+00 


O.OOOOOOOE+00 


P 


1 


-6.9546103E-02 


1.5249372E+00 


1.5944833E+00 


6.5786761E-01 


6.8791503E-01 


-1.9396281E-02 



This paper has been typeset from a TffjX/ WF^K file prepared 
by the author. 



